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Abstract 

The limited angular resolution of gamma-ray telescopes 
prevents the straight identification of the majority of the 
sources detected so far. While the only certified galac- 
tic gamma-ray emitters are associated with pulsars, more 
than 90% of the nearly 100 low latitude sources detected 
by EGRET on the GRO lack an identification. The best 
identification strategy devised over the years relies on the 
use of X-ray and optical data to single out possible coun- 
terparts, rating them on the basis of their extreme F x /F opt . 
Here we describe our multiwavelength programme based 
on the EPIC mapping of selected EGRET sources comple- 
mented by the optical coverage of the same fields obtained 
throug h the Wide Field Imager (WFI) of the ESQ 2.2 m 



gamma-ray error boxes (on average 1 sq deg) prevents a 
direct optical identification. For these reasons, the search 
for counterpart (s) has been pursued through the X-ray 
coverage of the fields followed by the optical study of each 
X-ray source. 

Figure 1 summarizes the "zoom in" approach devised dur- 
ing the 20 year long chase which led to the identification 



telesco pe . The field of view of the WFI, comparable to 
that of EPIC, will allow to directly cross-correlate X and 
optical data speeding up significantly the selection of in- 
teresting candidates worth follow-up optical (and X-ray) 
studies. 



of Geminga (Bignami & Caraveo 1996). The Einstein Ob- 
servatory was used to image the gamma-ray error box 
of Geminga, one of the brightest gamma-ray source in 
the galactic plane, shown in the upper right panel to- 
gether with the slighly fainter Crab pulsar. All but one 
of the four sources discovered by the Imaging Propor- 
tional Cou nter were readily optically identified ( Bignami 
ct al. 1983|) . The one lacking identification, 1E0630+178, 
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1. Introduction 



was later observed with the High Resolution Imager which 
yielded a few sq. arcsec error box. This was the target of 
CFH, Palomar and the ESO 3.6m telescopes, eventually 
yielding the my 25.5 G" as the candidate counterpart. 
This identification was later confirmed by the discovery 
of the source proper motion (Bignami ct al. 1993) which 
was found to improve the time solution of the gamma-ray 
photon (Mattox et al. 1996), providing a direct link be- 
tween the bright gamma-ray source and its faint optical 
counterpart. The link between the X and the gamma-ray 
behaviour was provided by the discovery of the pulsation 
in X-rays ( Halpern fc Holt 1992| ) followed by the confirma- 



The third EGRET catalogue ( |Hartman et al. 1999Q of 
high energy gamma-ray sources contains 271 entries as 
opposed to the 25 listed in the COS-B one. The low lat- 
itude (b < 10°), presumably galactic, sources went from 
22 to 80, while the high latitude (b > 10°), presumably 
extragalactic, ones jumped from 3 to 181, with blazars ac- 
counting for roughly half of them. 

Surprisingly enough, EGRET has done very little to clar- 
ify the nature of the low latitude, mostly galactic, gamma- 
ray emitters. While their number increased fourfold, suc- 
cessful identifications remain at a meagre < 10% level and 
rely only on the time signature of the gamma-ray pho- 
tons, which unambiguously associates these sources with 
pulsars. 

In the remaining cases, both the poorer statistics and the 
uncertain positions of the gamma-ray sources have ham- 
pered "blind" periodicity searches, while the size of the 



tion in the gamma-ray domain. The lack of radio emission 
prompted paraveo et al. (1996)| to classify Geming a as a 
radio quiet neutron star, the first of its kind, character- 
ized by a gamma-ray yield 1,000 times higher than the 
X-ray one, which , in turn, is 1,000 times higher than the 
optical one. Indeed, the identification procedure relies on 
the extreme F x /F opt value measured for Isolated Neutron 
Stars (INSs), 6 of which are kno wn to be gamma-ray emit- 
ters (see e.g. Thompson 2001 ). Such a multiwavelength 
approach is now the standard method to search for the 
counterparts of EGRET gamma-ray sources. 

Years of multiwavelength efforts have yielded (so far) 
less than a dozen tentative identifications encompassing 4 
energetic young radio pulsars, 3 Geminga-like radio quiet 
isolated neutron stars, 2 peculiar binary syst ems, 1 low lat- 



itude Blazar and few pulsars nebulae (see Caraveo 2001 
for a review of the proposed identifications). 



Proc. Symposium 'New Visions of the X-ray Universe in the XMM-Newton and Chandra Era' 26-30 November 2001, 



2 



Attempts to associate the low latitude unidentified EGRET 
sources with different classes of galactic objects (see e.g. 



applied to pinpoint the candidate neutron star countcr- 



part to the EGRET source 3EG J1835+5918 (Mirabal & 



Romeo et al. 1999; Gehrels et al. 2000) have not yielded, Halpern 2001) as well as to few other EGRET sources (see 



so far, gamma-ray source templates more appealing than 
the classical Isolated Neutron Star one. Thus, while not 
neglecting other possibilities, it is natural to assume that 
at least a fraction of the remaining low latitude sources 
are unidentified neutron stars. 

Our EPIC-ESO programme, centered on two middle lati- 
tude EGRET sources, is precisely aiming at these objects. 



|Caravco 2001 ). First of all, one has to start from a list of 



unidentified X-ray sources detected in the EGRET error 
boxes. Next step is to single out potential neutron star 
candidates taking advantage of their high F x /F opt as a 
distinctive character and using multicolor information as 
a further handle to solve ambiguous cases, e.g., when more 
optical entries are compatible with the X-ray position. 



2. INSS: RADIO-LOUD VS. RADIO-QUIET 

The search for neutron stars in gamma-rays started as 
soon as COS-B discovered the UGOs (Unidentified Gamma 
Objects, see Bignami & Hermsen 1983 for a review), but 
no new pulsar, after Crab and Vela, was unveiled. New 
searches have been spurred by the EGRET detections of 



3. Our XMM-ESO programme 

3.1. X-RAY SIDE 

Currently, we are focussing on two middle-latitude EGRET 
sources (3EG 0616-3310 and 3EG 1249-8330), selected on 
the basis of their relatively good positional accuracy, spec- 
tral shape, galactic location and lack of candidate extra- 



four more pulsars (e.g. [Thompson et al. 1996] ) but the lack galactic co unterpart. EPIC ( [Turner et al. 200l[ |5trudei 



substantially unchanged. Dedicated radio searches (Nice 



of results, experienced at the time of COS-B, appears to be |et al. 200l| ) is the ideal instrument to perform the X-ray 



fe Sayer 1997), aimed precisely at the search for radio pul- 



sars inside the error boxes of 10 of the brightest EGRET 
sources, yielded null results, showing that the straight- 
forward radio pulsar identification is not the only possi- 
ble solution to the enigma of the unidentified high-energy 
gamma-ray sources. This has been further strengthened by 
the work of |Nel et al. (1996) who investigated 350 known 



pulsars, finding few positional coincidences but no signifi- 
cant gamma-ray timing signature for any of the pulsars in 
the survey. The recently released portion of the Parkes sur- 
vey, encompassing 368 new pulsars, has been correlated to 
the EGRET error boxes, finding no more than 3 plausible 
pulsar candidates against a dozen of chance coincidences 
(Torres et al. 2001). With no gamma-ray instruments in 
operation, it is presently impossible to confirm these as- 
sociations. 

Indeed, gamma-ray astronomy does offer a remarkable ex- 
ample of an INS which behaves as a pulsar as far as X- 
and-gamma-astronomy are concerned but has little, if at 
all, radio emission. As an established representative of the 



non-radio- loud INSs (see Caraveo et al. 1996 for a review), 
Geminga offers an elusive template behaviour: prominent 
in high energy gamma-rays, easily detectable in X-rays 
and downright faint in optical, with sporadic or no radio 
emission. Although the energetics of a Geminga-like object 
is not adequate to account for the very low latitude (pre- 
sumably more distant) EGRET sources, the third EGRET 
catalogue contains several middle latitude sources which 
could belong to a local galactic population. In this case, 
their gamma-ray yield is certainly compatible with the 
rotational energy loss of a middle aged neutron star, like 
Geminga. 

Thus, it makes sense to apply to these sources the two- 
step strat egy devised during the 20 year long chase for 
Geminga ( Bignami fc Caraveo 1996 ), more recently also 



coverage of the gamma-ray error boxes, since it offers good 
angular resolution (< 6" FWHM), coupled to high sensi- 
tivity in a broad energy range (0.1-^12 keV) and good 
spectral resolution (E/AE w 20-j-50), over a wide field of 
view (sa 15 arcmin radius). Thus, each gamma-ray source 
error box, a circle of 30 arcmin radius, can be covered with 
four EPIC exposures. For a net exposure time of 10 000 
s for each pointing, this yields an homogeneous coverage 
of about 1 square degree down to a 5cr detection limiting 
flux ranging from w 7 x 10 -15 erg cm~ 2 s _1 to w 2 x 10 -14 
erg cm~ 2 s _1 , depending on the source spectrum. Taking 
these values as reference, folded with the LogN — LogS 



function for extragalactic sources (see e.g. Hasinger et al 
[2001 ) and including an evaluation of the galactic contri- 



bution (see e.g. puillout et al. 1996 ), we were expecting 
~ 100 sources (mainly AGN and active stars), positioned 
to within 5 arcsec, in each EGRET error box. To quickly 
identify such a number of serendipitous sources, we had 
to devise an ad hoc optical strategy. 

3.2. The optical side 

Given the range of F x /F opt values characteristic for the 
known classes of X-ray sources (see, e.g., Krauttcr et al 



1999), we ought to reach V ~ 25 in the optical follow-up in 



order to be able to discard most non-INSs identifications. 
Thus, although useful for a first filtering, the available 
Digital Sky Surveys are not deep enough for our purpose: 
dedicated optical observations are needed. 
However, the expected EPIC yield prompted us to seek an 
approach different from the one-by-one philosophy used so 
far to identify the few sources detected by previous X-ray 
telescopes in each gamma-ray error box. It is apparent 
that without a "massive" approach to the identification 
work the optical side is bound to became the bottleneck 
of our a multiwavelength chain. Our optical work should 
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be aimed at dozens of faint objects distributed over about 
1 sq deg. A dedicated, deep, multicolour optical survey of 
our fields must be the starting point for the real identifi- 
cation work. This is why we have decided to rely on the 
WFI of the 2.2 m as the optical complement of our EPIC 
coverage of EGRET sources. The WFI large field of view 
(33 x 34 arcmin), directly comparable to the EPIC one, 
will allow to cross correlate directly X and optical data 
speeding up dramatically the identification work. 

3.3. The state of the art 

The XMM/EPIC observation campaign of the two EGRET 
error boxes was completed in November 2001. The optical 
data are now being collected at the ESO 2.2m telescope 
(Period 68) . The analysis of the X-ray data is currently in 
progress, using the XMM-Newton Science Analysis Sys- 
tem (XMM-SAS) v5.2. 

Unfortunately, we discovered that space weather was not 
favourable during 2 of our 8 pointings. These observations 
(1 per EGRET error box) were badly affected by high par- 
ticle background episodes (the so-called soft proton flares, 
see XMM Users' Handbook) and we were forced to reject 
up to 80% of the integration time. 

We ran a preliminary source detection over the broad en- 
ergy band 0.3-8 keV, which gathers almost all the usable 
X-ray events from celestial sources. About 130 sources per 
EGRET error box were found and their fluxes, computed 
using a simple Geminga-like spectrum (a black body with 
kT=0.1 keV and a Nh of order 5xl0 20 cm -2 ), range from 



3 x 10 13 erg cm 2 s 1 



to 



10 14 erg cm 2 s 1 . 



Figure 2 shows the EPIC mosaic of the field of 3EG 1249- 
8330. The detected sources are marked with little circles. 
The South West observation clearly shows the dramatic ef- 
fect of flaring particle background. A more sophisticated 
source detection analysis is currently under way. Follow- 



ing Baldi et aL (2001), particular care is now devoted to 
create a background map which properly reproduce the 
background fluctuations on small spatial scales, a crucial 
step in order to assess the reality and the significance of 
our sources. 

As a first screening to filter out obvious non-INS iden- 
tifications, we take advantage of the available optical cata- 
logues to cross correlate our lists of sources (the position of 
which is known within 5 arcsec). Our starting point is the 
extended version of the recently released Guide Star Cat- 
alogue 2 (GSC2), an all-sky, multi-epoch and multicolor 
optical catalogue based on photographic surveys carried 
out between 1953 and 1991 (McLean et al. 2002-in prepa- 
ration). The GSC2 provides color information in at least 
three photographic passbands (roughly corresponding to 
B, R and /) down to B ~ 22 and morphological clas- 
sification at a ~ 90% confidence level for objects in our 
latitude range and brighter than B ~ 19. To extend our 
color coverage to the infrared (J, H and K) we use the 
updated release of the 2 Micron All Sky Survey (2MASS) 



catalogue. 

After the cross correlations, the color catalogues are com- 
bined and used as input for the object classification proce- 
dure, with the morphological classification from the GSC2 
as a further aid. For this, we apply the automatic algo- 



rithm developed by Hatziminaoglou, et al. 200C , which fits 



model fluxes, simulated from template spectral libraries, 
to the observed ones through a % 2 minimization technique. 
Th e reliability of this method h as been tested successfully 
by Hatziminaoglou, et al. 2002] for the classification of the 
objects detected in the Chandra Deep Field South. Fi- 
nally, the properties of the X-ray sources (e.g. hardness 
ratios, spectra) will be compared with the classifications 
of their candidate optical counterparts and identifications 
evaluated through a decision tree. 

The whole procedure will be then repeated for the remain- 
ing, non identified, EPIC sources, using deeper optical cat- 
alogues extracted from the WFI data. The observations 
(in UBVRI) are now being carried out in Service Mode 
by the 2.2m team. So far, only the field of 3EG 0616-3310 
has been observed, with an area coverage of ~ 25%. 
Of course, only the availability of the multicolor optical 
data from WFI will allow us to perform the final step of 
our identification chain. 



4. Conclusions 

The identification of gamma-ray sources, both individu- 
ally and as a population, rests on multiwavelength obser- 
vations. The programme we are pursuing at X-ray and op- 
tical wavelengths will yield, in a relatively short time, can- 
didate counterparts worth follow-up investigations. More- 
over, we plan to make available the complete catalog of 
the sources detected in our 1.6 sq. deg medium galactic 
latitude survey. 
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Figure 1. From upper right to lower left: sketch of the 20 years-long identification chain which led to the identification of 




Figure 2. EPIC mapping of the error box of 3EG 1249-8330. North is up, East is left. The detected sources, marked 
with circles, are more than 130. The lower right pointing is the one affected by flaring background. The reduction of 
the fraction of good observing time is evident. 



